arXiv:1506.03386v2 [hep-ph] 23 Sep 2015 


TUM-HEP 1000/15 


A novel approach to derive halo-independent 
limits on dark matter properties 


Francesc Ferrer 1 , Alejandro Ibarra 2 , Sebastian Wild 2 

1 Physics Department and McDonnell Center for the Space Sciences, 
Washington University, St Louis, MO 63130, U.S.A. 

2 Physik-Department T30d, Technische Universitat Miinchen, 
James-Franck-Strafie, 85748 Garching, Germany 

September 24, 2015 


Abstract 

We propose a method that allows to place an upper limit on the dark matter elas¬ 
tic scattering cross section with nucleons which is independent of the velocity distribu¬ 
tion. Our approach combines null results from direct detection experiments with indi¬ 
rect searches at neutrino telescopes, and goes beyond previous attempts to remove as- 
trophysical uncertainties in that it directly constrains the particle physics properties of 
the dark matter. The resulting halo-independent upper limits on the scattering cross 
section of dark matter are remarkably strong and reach cigj < 1CP 43 (10 -42 ) cm 2 and 
<7gD 1$ 10 -37 (3 x 10 -37 ) cm 2 , for dark matter particles of tudm ~ 1 TeV annihilating into 
W + W~ (bb), assuming p\ oc = 0.3 GeV/cm 3 . 


1 Introduction 

Save for its gravitational effects, which establish dark matter (DM) as the second largest con¬ 
tributor to the energy budget of the Universe [I], very little is known about its composition 2 jb]. 
An attractive candidate for making up the dark matter is a weakly interacting massive particle 
(WIMP), that could be unveiled through its weak-scale interactions. 

Two complementary approaches have been proposed to probe the hypothetical WIMP pop¬ 
ulation inside the Solar System. Direct detection experiments aim to detect the nuclear recoil 
induced by the elastic scattering of the dark matter particles traversing a detector at the 
Earth |6|. The second approach consists in the search using neutrino telescopes of the high 
energy neutrinos which are hypothetically produced in the annihilation of dark matter parti¬ 
cles which have been previously captured in the Sun via a series of scatterings with the solar 
matter |7|. 

Over the past few years, these programs have attained the sensitivity to probe the particle 
physics properties in a vast range of WIMP scenarios. The interpretation of these experimental 
efforts, however, is hampered by uncertainties in the required astrophysical input. For instance, 
the event rates at both neutrino telescopes (NT) and direct detection (DD) experiments are 
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sensitive to the elastic DM-nucleon scattering cross section. But to translate experimental 
data into upper limits on the DM scattering cross section requires the input of a DM velocity 
distribution f(v). 

A common strategy is to adopt a Maxwell-Boltzmann velocity distribution in the galactic 
rest frame. Using this scheme, upper limits on the spin-independent and spin-dependent WIMP- 
nucleon scattering cross section have been published by direct detection experiments, such as 
LUX §, XENONIOO §, SuperCDMS |10|, PICO (IT), COUPP [12] and SIMPLE [l3j, as 
well as by neutrino telescopes, such as Super-Karniokande (l4]| and IceCube 15 . A Maxwell- 


Boltzmann velocity distribution, however, only holds when the dark matter distribution is 
modeled after the Standard Halo Model (SHM) [16], an isothermal sphere profile with an 
isotropic velocity distribution. Recent N-body simulations, however, provide evidence that 
a Maxwellian distribution is not a good description of the smooth halo component 17-19 


Furthermore, the dark matter halo of our Galaxy might contain tidal streams or a dark disk 
component 20 -23 , adding a new source of uncertainty to the dark matter velocity distribution 
inside the Solar System. Therefore, to robustly constrain the WIMP scattering cross section 
from experiments, it is necessary to devise methods that do not depend on the choice of the 
velocity distribution. 

Several methods have been proposed to derive halo-independent statements from direct 
detection experiments 24-36 . Broadly speaking, the rationale behind these studies can be di¬ 


vided into two categories. First, the experimental data can be expressed in terms of a variable 
that includes both the information about the DM scattering cross section and the (integrated) 
velocity distribution, often denoted as By comparing measurements and upper lim¬ 

its of such an observable from different direct detection experiments, one can quantify the 
compatibility of a positive claim with a null result in a halo-independent way. Second, the 
existing methods can be used to infer fundamental particle physics properties of DM from pos¬ 
itive signals in current or future experiments. Along these lines, it has also been shown that 
the complementarity of direct searches and neutrino telescopes can be important, due to the 
different dependence of both detection methods on the DM velocity distribution 


37,38 


In this work we propose a novel method which allows for the first time to use existing 
data for calculating an upper limit on the scattering cross section which is independent of the 
velocity distribution, by combining null results from direct detection experiments and neutrino 
telescopes. 

The paper is organized as follows. In section [2] we outline our method and in section [3] we 
derive an upper limit on the elastic scattering cross section from a direct detection experiment 
or from a neutrino telescope for a general velocity distribution. Then, in section [4] we apply 
our method to derive an upper bound on the cross section which is independent of the velocity 
distribution. This is the main result of this paper. Lastly, in section [5j we apply our method to 
derive a halo-independent lower limit on the scattering cross section from a hypothetical signal 
in a future direct detection experiment. We present our conclusions in section [6j and provide 
details of direct detection data in appendix [A] 


2 Probing WIMP dark matter inside the Solar System 

We assume that the Solar system is embedded in a distribution of WIMP dark matter particles 
with mass density p\ oc . The DM halo is spatially homogeneous on Solar system scales and its 
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velocity distribution f(v) in the rest frame of the Sun is normalized as 


d 3 vf{v) = 1 , 


( 1 ) 


>V<V n 


where v = |u| and u max is the maximal velocity of the dark matter particles in the galactic halo. 
Under the common assumption that all dark matter particles in the halo are gravitationally 
bound to the Galaxy, u max ~ 777 krn/s, which is the sum of the galactic escape velocity ~ 533 
krn/s ||39j and the local velocity of the Sun with respect to the halo ~ 244 km/s (40 421. We 
comment in section [4] on the implications of relaxing this assumption. 

The number of expected recoil events at a DD experiment can be expressed as: 


R = £-Y' [°° dE R e(E R ) iiPloc [ 

'o m Ai m DM J v>v vm (E) 
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Here, v denotes the dark matter velocity in the detector frame, hence the velocity distribution of 
dark matter particles is f(v + v 0 bs(t)), with v ohs (t) = u©(f) being the velocity of the Earth with 
respect to the solar frame, and |u©(f)| = 29.8 km/s (we note again that we define f{v) in the rest 
frame of the Sun). Furthermore, dai/dE R is the differential scattering cross section of a WIMP 
off a nuclear isotope i with mass m Ai and mass fraction G in the detector, and {E R ) = 


^Jm Ai E R / (2//^.) is the minimal speed necessary for a dark matter particle to induce a recoil 
with energy E Rj with ji Ai being the reduced mass of the WIMP-nucleus scattering. Lastly, 
e(E R ) and £ are the detection efficiency and exposure, respectively, both depending on the 
specific experiment under consideration. 

On the other hand, the neutrino fiux from annihilations in the Sun is completely determined 
by the capture rate, assuming that dark matter capture and annihilation are in equilibrium. 


The capture rate is given by 43 
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where r]i(r) is the number density of the element i at a distance r from the solar center, u eS c(^) 
is the escape velocity, and = 2 v esc (r) — rri Al is the maximum speed 

of a dark matter particle such that the capture in the Sun remains kinematically possible. 
For scattering cross sections that can be currently probed by neutrino telescopes, equilibration 
is achieved as long as the annihilation cross section multiplied by the relative dark matter 
velocity satisfies (cm) ann > 1CT 28 cm 3 /s In the rest of this work, we will always assume that 
equilibrium has been reached. 

The largest uncertainties in the calculation of the scattering and the capture rates in a 
given particle physics model stem from our ignorance of the WIMP scattering cross section 

1 Capture and annihilation in the Sun are in equilibrium if a eq = tanh 2 ( tQ-\/C ■ 14) is close to one, where 
f Q ~ 1.5 • 10 17 s is the age of the Sun, and 14 ~ 5.17 • 10 -57 t • 3 . 10 i 26 J c m 3 / s ' (lleV annihilation 

constant [I]. It is straightforward to check that for (av) ann > 10 -28 cm 3 /s and for the capture rates C that are 
currently probed by IceCube or Super-Kamiokande, one has a e q <1 0.99. 
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with nuclei and of the dark matter velocity distribution. It is common in the literature to cast 


the differential cross section as 44 


da* 
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where ctsi and usd are, respectively, the spin-independent and spin-dependent cross sections at 
zero momentum transfer, which can be calculated in a concrete dark matter model in terms of 
its fundamental parameters, while F h c,-[(E ri ) and are form factors that depend on the 

nucleus, and which are given in [45 and in [46] for the spin-independent and the spin-dependent 
interactions, respectively. Furthermore, most analyses adopt a Maxwell-Boltzmann velocity 
distribution in the galactic rest frame. Under this assumption, the null searches for dark matter 
particles in the Solar System using direct detection experiments and from neutrino telescopes 
are translated into limits on crsi and Usd, and confronted with the predictions from models. 
Clearly, a particle physics model violating these limits is not necessarily excluded, since the 
actual velocity distribution might deviate from the simple Maxwellian form. It is then important 
to derive limits on asi and Osd which do not depend on the form of the velocity distribution, in 
order to robustly exclude a concrete particle physics scenario. In the following sections we will 
present a novel method to obtain halo-independent limits on the cross section which exploits 
the complementarity between direct detection experiments and neutrino telescopes. 


3 Dark matter in the solar system as a superposition of 
streams 

For the purpose of deriving limits on the WIMP cross section for general velocity distributions 
we find it convenient to decompose the dark matter population in the solar system as 

f(v)=l d 3 v 0 d {3) {v-v 0 )f(v 0 ), (5) 

J |tTo|<v max 

which physically can be interpreted as a superposition of hypothetical streams with fixed ve¬ 
locity v 0 with respect to the solar frame, fy 0 (v) = S^(v — v 0 ). Then, using Eqs. ([ 2 J and (J 3 ]) , 
the number of expected scattering events R and the capture rate C can be cast as 

R = I d 3 U 0 / (vo) Rvo , 

J |ub | <-« max 

C= I d 3 u 0 /(F 0 )CV 0 , (6) 

J\vo\<v max 

where R$ 0 and C$ 0 are, respectively, the number of scattering events and the capture rate for 
the dark matter stream with velocity distribution fv 0 {v ) = d^ 3 \v — Vq ). 

For every stream with velocity vq with respect to the Sun (and therefore with velocity 
Vq — Ve(R) with respect to the Earth) one can obtain an upper limit on the scattering cross 
section from the null results of a direct detection experiment by requiring that R$ 0 < i? max - 
Here, i? max is the upper limit on the expected number of scattering events following from the 
results of the corresponding experiment. For definiteness, we assume equal coupling to protons 
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and neutrons for the spin-independent scattering and coupling only to protons for the spin- 
dependent scattering. Appendix [A] provides further details of the calculation. We show in 
Fig.0 for illustration, the 90% C.L. upper limits on the spin-independent (upper panels) and 
spin-dependent cross sections (lower panels) for tudm = 100 GeV as a function of the stream 
speed Vq = % which follow from the null results from XENONIOO or COUPP, respectively. 
The dashed red lines show the limits for various angles between the stream velocity % and 
a fixed Earth velocity ve, while the red solid line corresponds, for a fixed speed vq, to the 
weakest among the limits for all angles. This conservative upper limit on the cross section will 
be denoted as (Xma^A)), either for the spin-independent or the spin-dependent cross section. 
Then, by construction, Rff 0 (a) > _R max for a > cr ^ (r> 0 ). 

Similarly, for every stream there is an upper limit on the scattering cross section which is 
allowed by a neutrino telescope, which we calculate from the requirement that CV 0 < C max , 
where C ma , x is the maximum capture rate allowed for a given annihilation channel, assuming 
equilibration between capture and annihilation. For the calculation of the capture rate induced 
by the spin-independent coupling we include the 29 most abundant elements in the Sun and 
we assume Gaussian form factors, as parametrized in 47 . On the other hand, to obtain 


the capture rate induced by the spin-dependent coupling to the proton we take into account, 
not only the scattering off hydrogen, as commonly done in the literature, but 

For the number density profile we adopt the 


following 48 


48 


on 


also off 14 N, with form factors as given in 

solar model AGSS09 |49|. The blue solid line in Fig. [I] shows the 90% C.L. upper limits 
the spin-independent (upper panels) and spin-dependent cross sections (lower panels) for dark 
matter streams with speed Vq = 1% which follow from the requirement that the capture rate 
in the Sun is smaller than the maximum value C max inferred from the IceCube data. We take 
again uidm = 100 GeV, and consider the annihilation final states W + W~ (left panels) and bb 
(right panel). The noticeable kink in the IceCube limit on the spin-dependent cross section is 
due to the fact that at low velocities the capture rate is dominated by scatterings off H, while 
scatterings off 14 N become important at large velocities. In the case of neutrino telescopes, the 
most conservative upper limit on the cross section for a given stream speed Vo which follows from 
the null results of a neutrino telescope will be denoted as cr^J x (n 0 ), therefore, C%(a) > C max 
for a > cr£J x (w 0 )- 

The limits on the scattering cross section from direct detection or from neutrino telescopes 
derived for streams with speed t’o can be used to calculate limits for a general velocity distri¬ 
bution f(v). Noting that the differential cross section is linear in ctsi, ct sd , it follows that the 
number of expected scattering events satisfies 


Rv 0 ( a ) DD ( \-^h)[ a max(' l; o)] A DD / \ 

^maxWW u maxv^ 


-max 


(7) 


where a denotes either crsi or <7 sd, and we have used the fact that, for the dark matter stream 
with velocity v 0 , cr^ x (n 0 ) gives the largest scattering rate allowed by the direct detection 
experiment. Inserting this inequality in Eq. ([6]) gives 
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Finally, requiring that the number of expected events induced by the cross section cr is in 
agreement with the experimental upper limit, i?(cr) < i? max , one obtains from a given direct 
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Peak position Uq of the velocity distribution 


Upper limit on tr^, for fv 0 (v) = 5 3 (v — Uq) 





Peak position vq of the velocity distribution Peak position vq of the velocity distribution 


Figure 1: Upper limits on the spin-independent dark matter-proton scattering cross section <7gj, nor¬ 
malized to a local dark matter density p\ oc = 0.3 GeV/cm 3 , derived from the null results of XENONIOO 
and IceCube (upper plots) and on the spin-dependent dark matter-proton scattering cross section <jg D 
from the null results of COUPP and IceCube (lower plots), assuming annihilations into W + W~ (left 
plots) or bb (right plots), for todm = 100 GeV and for a velocity distribution corresponding to a stream 
with speed vq with respect to the rest frame of the Sun, fv Q (v) = 5^ 3 \v — vq). The various dashed 
red lines denote the upper limits for different angles between vq and the velocity of the Earth. In 
each plot, the maximal value of the cross section allowed by the direct detection experiment and the 
neutrino telescope is denoted by <r*. 
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detection experiment the upper limit on the cross section 
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A similar rationale applied to the capture rate at a neutrino telescope gives 


a < 


-J\vo\<v 
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d 3 f(Vo) 

'C(«<0 


n -1 


(9) 


( 10 ) 


This formalism allows to readily calculate upper limits on the cross section for a general ve¬ 
locity distribution from the function cr max (v 0 ) derived for a given direct detection experiment 
or neutrino telescope, without the necessity of going through the rather involved calculations 
described in Section 00 


4 A halo-independent upper limit on the scattering cross 
section 

Direct detection experiments and neutrino telescopes probe the WIMP population in the Solar 
System in a complementary way: direct detection experiments are insensitive to slow moving 
WIMPs, since the recoiling nucleus has an energy below the detector threshold. However, 
this population of dark matter particles can be efficiently captured inside the Sun and thus 
leads to a high energy neutrino flux which could be detected in a neutrino telescope. This 
complementarity is apparent from Fig. [l]for dark matter streams, since for every stream speed 
Vo there exists a finite upper limit on the cross section, which follows either from a direct 
detection experiment or from a neutrino telescope. 

We will show in this section that the complementarity between direct detection experiments 
and neutrino telescopes makes it possible to derive an upper limit on the scattering cross which 
is independent of the velocity distribution by using the upper limits on the cross section for 
dark matter streams. To this end, we first determine, for a given dark matter mass and a given 
annihilation channel, the largest value of the scattering cross section cr* which is allowed by a 
concrete direct detection experiment and a concrete neutrino telescope, assuming a stream-like 
velocity distribution with speed between 0 and t> max : 

a, = max (u), (u max )} , (11) 

where v is the speed for which cr aiax (u) = cr ^ Y (u). The determination of the maximal scattering 
cross section allowed by a direct detection experiment and by a neutrino telescopes is illustrated 
in Fig. [l] for a dark matter mass m^M = 100 GeV which annihilates into W + W~ (left panels) 
and bb (right panels), considering XENONIOO and IceCube for the spin-independent limits 
(upper panels), and IceCube and COUPP for the spin-dependent limits (lower panels). Then, 
by construction, 


2 The functions cr laax 


O-maxK) < a * for V < W < t'max , 

o^ivo) < a* for 0 < v 0 < v . 

(vo) and crSaxf^o) are available from the authors upon request. 


( 12 ) 

(13) 
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To derive an upper limit on the cross section by combining results from a direct detection 
experiment and a neutrino telescope we first note from Eq. Q that 
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Using now Eq. (12) and defining 


Sf = 


d 3 u 0 f{v 0 ) , 


we obtain 




< -r 
°f 


(14) 


(15) 


(16) 


An analogous calculation applied to the upper limit on the cross section from a neutrino 


telescope, Eq. (10), gives, using Eq. (13), 


a ^ T~T“ » ( 17 ) 

1-df 

where the factor (1 — 5/) follows from the normalization of the velocity distribution, Eq. (jlj), 
and the definition of 5/, Eq. (15). 

Finally, the upper limits on the cross section from a direct detection experiment and from 
a neutrino telescope, Eqs. (16) and (0, imply 

a < 2<j* . (18) 

This is the main result of this paper: for a given dark matter mass and a given annihilation 
channel, there exists an upper limit on the scattering cross section which follows from com¬ 
bining the null results from direct detection experiments and neutrino telescopes and which 
is independent of the choice of the halo velocity distribution f(v). In this equation, cr* can 
be calculated following the procedure described in Section [2] from the analysis of dark matter 
streams. It can also be checked that the distributions that saturate the limit Eq. (18) satisfy 
5 , = 1/2, 

The halo-independent upper limits on the spin independent and spin dependent cross sec¬ 
tions are shown in the upper and lower plots of Fig. [2j respectively, as a function of the dark 
matter mass, normalized to a local dark matter density p\ oc = 0.3 GeV/cm 3 . The left panels 
correspond to dark matter annihilating into W + W~ /t + t~ , while the right panels assume an¬ 
nihilation into bb. For each dark matter mass we calculate, using the method described above, 
an upper limit on the cross section using the null results from one direct detection experi¬ 
ment (either XENONIOO or SuperCDMS for the spin-independent scattering, or COUPP or 
SIMPLE for the spin-dependent scattering) and one neutrino telescope (either IceCube or Su- 
perKamiokande). Finally, we take the upper limit on the spin-independent and spin-dependent 
cross sections which is allowed by all the experiments considered in our analysis. 


3 We do not consider limits from LUX [§], as the information provided by the collaboration is not sufficient 


to derive limits for non-standard velocity distributions 
experiment 


50 . Similarly, also the recent results of the PICO 


11 can not be generalized to other velocity distributions without additional information about 


efficiency of the detector. 
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Figure 2: Halo-independent upper limits on cjgj (upper panel) and 0 g D (lower panel), normalized to 
a local dark matter density p\ oc = 0.3 GeV/cm 3 , for annihilation into W + W~ /t + t~ (left panels) and 
bb (right panels). The red solid curves correspond to v max = 777 krn/s, while the red dotted curves 
assume v max = 0.05 c. Also, the black solid and dashed curves show the upper limits from different 
direct detection experiments and neutrino telescopes, assuming the Standard Halo Model (SHM) as 
defined in the text. 

In order to compare our results with the standard exclusion bounds from direct detection 
experiments as well as from neutrino telescopes, we also show in Fig. [2] the upper limits on the 
scattering cross section for the specific case of the Standard Halo Model (SHM), i.e. assuming 
a Maxwell-Boltzmann velocity distribution with velocity dispersion u 0 — 230 km/s and galactic 
escape velocity v esc = 533 km/s. In the upper panels, corresponding to spin-independent 
scattering, black solid lines show the limits from Super-CDMS and XENONIOO, while in the 
lower panels, corresponding to the case of spin-dependent scattering off protons, they show 
bounds from SIMPLE and COUPP. Besides, the black dashed lines show the upper limits from 
Super-K and IceCube. The halo-independent limits are, as expected, somewhat weaker than the 
combined limit from experiments assuming the standard halo model. Nevertheless, our limits 
are remarkably strong and reach Ugj < 1CT 43 (10~ 42 ) cm 2 and (Jg D < 1CP 3 ' (3 x 1CF 37 ) cm 2 , for 
annihilations into W + W~ (bb) at tudm = 1 TeV, assuming p\ oc = 0.3 GeV/cm 3 . For the spin 
independent coupling to protons these limits are better than those obtained by IceCube for the 
SHM, while for the spin dependent coupling to protons, only a factor of a few worse than the 
limits from the SIMPLE or COUPP experiments, also assuming the SHM. 
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The limits presented in this work rely on rather weak assumptions. First, we assume that 
the dark matter density and velocity distribution at the position of the Sun and the Earth 
are identical and constant over the equilibration time of captures and annihilations, which we 
estimate to be ~10-100 million years or about 0.2-2% of the age of the Sun. We also assume 
Umax = 777 km/s, which follows from the well motivated assumption that all dark matter 
particles in the Solar System are gravitationally bound to the galaxy. Nevertheless, we have 
checked that even if the value of the maximal speed is doubled, our limits do not change 
significantly, since the largest value of the scattering cross-section er* is typically attained at 
the velocity v, unless u max is very large, as can be seen from Fig. [I] (we recall that v is defined 
as the velocity for which o~ ^ aY (n) = cr ^J Y (h)). To include the possibility of a very large maximal 
velocity, we have also calculated the corresponding limits for v max = 0.05c, which approximately 
corresponds to the largest value of the velocity where the standard non-relativistic description 
of the scattering and capture rates in Eqs. i§. i§ can be applied. The limits are shown in 
rig- i as a dotted line and are, at most, one order of magnitude worse than those derived for 
tVa = 777 km/s. 

The halo-independent upper limits on the scattering cross section presented in this work 
will improve with upcoming data from direct detection experiments and neutrino telescopes, 
assuming no signal is detected. It follows from Fig. [l]that the limits on the spin-independent 
cross section can be improved by decreasing the threshold of the corresponding direct detection 
experiment, or by increasing the sensitivity of neutrino telescopes. In the near future, this 
could be achievable with e.g. Super-CDMS SNOLAB 51 and the results from the 86-string 
configuration of IceCube [52]. On the other hand, for spin-dependent interactions, the current 


limiting factor is the exposure of the relevant direct detection experiments, as can be seen from 
the lower panels in Fig. [I] Here, future upgrades of the PICO experiment 
improving the halo-independent upper limits. 


53 would help in 


5 Lower limit on the cross section from signal events 


We briefly discuss in this section the implications for the determination of the dark matter pa¬ 
rameters in the case of a positive signal in a xenon based experiment using our halo-independent 
approach. In contrast to the previous discussion, in this section we do not consider the com¬ 
plementary information arising from neutrino telescopes. We consider for concreteness an 
experiment with an exposure of 1 ton ■ yr, sensitive to the energy range between 3 and 45 keV, 
with 100% detection efficiency and an energy resolution cr(E R ) = 0.6 keV \/E R /keV 54 . We 
also consider a benchmark scenario where the true dark matter parameters are m{^ c ^ = 100 


GeV and <Tg[ tluc) = 10 -46 cm 2 , which approximately lies one order of magnitude below the cur¬ 
rent LUX limit. Assuming the SHM as the true velocity distribution, this benchmark scenario 
would lead to ~ 19 events in the future experiment, which amounts to a 90% C.L. lower limit 
on the number of recoil events R min = 14.0. 

On the basis of this number of events, we now calculate a halo-independent lower limit on 
the scattering cross section following a similar rationale as in Section [3] We first determine, for 
a given dark matter mass, a lower limit on the scattering cross section assuming a dark matter 
stream with velocity v 0 with respect to the Sun from requiring R$ 0 > R min , and we denote 


VO 

a™(v 0 ) the strongest among the limits for all possible angles between the stream velocity v 0 


and a fixed Earth velocity v R . Then, by construction, R^ q {<j ) < R m in for cr < cr°°(u 0 ). This 
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Figure 3: Lower limit on <7gj, normalized to a local dark matter density p\ oc = 0.3 GeV/cm 3 , derived 
from the hypothetical detection of 19 recoil events at a xenon based direct detection experiment, for 
tttdm = 300 GeV and for a velocity distribution corresponding to a stream with speed vq with respect 
to the rest frame of the Sun, fv 0 (v) = 5^ 3 \v — vq ). The minimum value of the cross section leading to 
the hypothetical number of recoil events is denoted by a'. 


procedure is illustrated in Fig. [3] for tudm = 300 GeV. The dashed red lines show the limits 
for various angles between the stream velocity vq and a fixed Earth velocity ve, while the red 
solid line corresponds to cr°P(u 0 ). Finally we construct, following similar steps as in Section 
[3j a lower limit on the scattering cross section for a general velocity distribution from the 
corresponding limits for the stream-like distributions: 

(19) 


a > 


d 3 


u 0 - 


/(To) 


VlwiKu 1 


a, 


DD 


Oo). 


For each dark matter mass the function cV|P (u 0 ) has a minimum, which we denote as a' (see 
Fig. g. Then, from Eqs. (|TJ) and ( |l9j ) we obtain the following lower limit on the scattering 
cross section, which is independent of the velocity distribution: 


a > a' . 


( 20 ) 


Fig. [4] shows, as a red solid line, the lower limit on the spin independent scattering cross 
section with protons as a function of the dark matter mass, calculated using this method. For 
comparison, the Figure also shows, as a gray band, the reconstructed 90% C.L. region assuming 
the SHM, as well as the halo-independent lower bound recently proposed in 55 based on the 
general inequality 


J U>«mm 


1 

5 

^min 


( 21 ) 


where u m i n is the minimal speed necessary for a dark matter particle to induce a recoil with 
energy above the threshold of the experiment (c/. Section [2]). Notably, for the true dark matter 
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Figure 4: Halo-independent lower limit on cjgj, normalized to a local dark matter density p\ oc = 
0.3 GeV/crn 3 , from the hypothetical observation of 19 recoil events in a xenon based direct detection 
experiment (red line), reconstructed 90% C.L. region assuming the SHM (gray band), and halo- 
independent lower limit derived from Eq. (21) (blue dotted line). The black dot indicates the true 
dark matter mass and scattering cross section. 


mass, the lower limit lies only a factor of ~ 2 below the true cross section (shown as a black 
point). We have checked that this conclusion does not change significantly with the number of 
observed events, as long as it is larger than ~ 5. Furthermore, had the true dark matter mass 
been different, for the same number of signal events one finds a lower limit on the cross section 
which lies only a factor ~ 2 — 3 below the true value, provided > 30 GeV. This is due 

to the degeneracy between the dark matter mass and the scattering cross section in the event 
rate, which is reflected for the SHM in the gray band in Fig. [4j 

Finally, we note that our method only uses the information about the total number of 
events, and hence it is particularly well-suited when only a few events are detected. On the 
other hand, once enough data is available to infer the recoil energy spectrum, halo-independent 
methods using the spectral information can break the degeneracy between the dark matter 


mass and scattering cross section 33,56 


6 Conclusions 

We have presented a new method to calculate limits on the dark matter scattering cross section 
for a general dark matter velocity distribution. Our method is based on decomposing the 
velocity distribution in dark matter streams with fixed velocity. We have then calculated the 
upper limit on the scattering cross section as a function of the stream speed for various direct 
detection experiments and neutrino telescopes and we have presented a simple formula that 
allows to calculate the corresponding limit for a general velocity distribution, based on this 
decomposition. This formula can be useful for a fast calculation of the limits for a general 
velocity distribution, since all the details of the calculation of the number of recoil events for a 
given direct detection experiment or of the capture rate in the Sun are already included in the 
limits for the stream distributions. 
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We have applied this method to calculate an upper limit on the dark matter scattering 
cross section by combining the null results from direct detection experiments and from neu¬ 
trino telescopes and which is independent of the dark matter velocity distribution. Our only 
assumptions are that the dark matter density and velocity distribution at the position of the 
Sun and the Earth are identical and constant over the last 10-100 million years, and that the 
dark matter particles move non-relativistically. The resulting limits are remarkably strong and 
reach ergj < 10 -43 (10 -42 ) cm 2 and <jg D < 10 -37 (3 x 10~ 37 ) cm 2 , for annihilations into W + W~ 
(■ bb ) at ??tdm — 1 TeV, assuming p\ oc = 0.3 GeV/cm 3 . 

Lastly, we have presented a simple procedure to derive a halo-independent lower limit on 
the scattering cross section as a function of the dark matter mass in the case of the detection of 
signal events in a future experiment. We have illustrated this procedure for a ton-scale xenon 
based experiment and we have found a lower limit which lies, for a true dark matter mass larger 
than ~ 30 GeV, only a factor of a few below the true cross section. 
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A Calculating upper limits from direct detection data 


In this appendix we provide the details of our method of deriving upper limits on the scattering 
cross section from XENON100, SuperCDMS, COUPP, and SIMPLE. For each of these experi¬ 
ments, we obtain the upper limit by comparing the number of expected recoil events R , given 
by Eq. (J2]), with the maximal number of events i? max allowed by the result of the experiment. 

More precisely, 


For XENON100, we calculate the detection efficiency e(E R ) following 57 


we use the scintillation efficiency C e g(E R ) from 58 , conservatively setting it to zero below 


E r = 3 keV. The detection efficiency as a function of the SI signal, as well as the efficiency 
related to the S2 threshold cut is taken from |9j. After running with an exposure of 7636.4 
kg • days, XENON100 observed two events in the predefined signal region given by S'™ 111 = 3 PE 
and S™' dX = 30 PE. Making no assumption about the background, the 90% C.L. upper limit on 
the number of recoil events is given by i? max = 5.32. 

in the recoil 


The SuperCDMS limits are obtained by using the efficiency e{E R ) given in 10 


energy range between 1.6 keV and 10 keV. With an exposure of 577 kg-days, SuperCDMS 
observed 11 candidate events. As there is no full understanding of the expected background, 
we derive limits by conservatively assuming zero background events, leading to if max = 16.6 at 
90% C.L. 


The latest result of the COUPP collaboration [12 is presented in the form of three non- 
overlapping data sets with different nucleation thresholds, given by 7.8 keV, 11.0 keV, and 
15.5 keV. The effective exposures (number of observed events) are 55.8 kg-days (2 events), 
70.0 kg-days (3 events), and 311.7 kg-days (8 events), respectively. Making no assumption 
about the background, the maximal number of expected events at 90% C.L. is given by _R max = 
5.32, 6.68, and 12.99, respectively. We derive individual upper limits for each of the three data 
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sets, choosing the most constraining one separately for each dark matter mass and interaction 
type. In all cases, the efficiency is assumed to be 0.49 above threshold for scattering off C or 
F, and 1 for scattering off I jl21; furthermore, we only consider recoils up to 500 keV. 

Finally, we derive limits from the SIMPLE experiment by using the combined Stage I and 
II results presented in 13 . Following the method pursued by the collaboration, we use an 


energy threshold of Et = 8 keV, with an efficiency e(En) = 1 — exp (—3.6 (Er/Et — 1)). As 
for COUPP, we limit ourselves to recoil energies below 500 keV. After running with an effective 
exposure of 20.18 kg-days, 11 events were observed. With the same conservative assumption 
b = 14.53 for the expected number of background events as used in 


13 


one obtains f?, max = 4.02 


at 90% C.L, using the Felclman-Cousins procedure 59 


For all experiments considered in this work, we checked that, under the assumptions specified 
above, the upper limit derived using the standard Maxwell-Boltzmann distribution matches 
sufficiently well the limit published by the corresponding collaboration. 
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